Unique Translation between Hamiltonian Operators and Functional Integrals 
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A careful treatment of the discretization errors in the path integral formulation of quantum 
mechanics leads to a unique prescription for the translation from the Hamiltonian to the action in 
the functional integral. An example is given by an interaction quadratic in the occupation number, 
characteristic for manybody bosonic systems. As a result, the term linear in the occupation number 
(chemical potential) receives a correction as compared to the usual formulation based on coherent 
states. A perturbative calculation supports the relevance of this correction. 
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I. INTRODUCTION 

Many modern techniques in particle and statistical 
physics rely on theories given in the form of functional in- 
tegrals, while the basic description of the system is often 
easier and more intuitively formulated on the operator 
level. The transition between these two languages has 
been given much attention (see e.g. 0-^) and is recog- 
nized as a mathematically complicated issue if the Hamil- 
tonian cannot be separated into a momentum and a lo- 
cation dependent part. Nevertheless, calculations based 
on functional integral formulations have been very suc- 
cessful in applications, for instance in the context of the 
renormalization group in statistical physics. 

In order to interpret the parameters which appear in 
the functional integral correctly, it is necessary to con- 
nect them to the parameters of the Hamiltonian opera- 
tor Ti. This translation between the parameters in the 
functional integral and the Hamiltonian has to be unique 
since physical quantities can be computed both ways. (If 
necessary, manybody models have to be regularized such 
that the expectation values of Ti and other operators are 
well defined. Typically, this can be done by discretization 
on a lattice.) 

For an ensemble of interacting bosons of mass m with 
the Hamiltonian 



.,v{q3)aq^aq^, (1) 



?i,92,g3 



the partition function is written as a functional integral 
over a complex field ^Pq{T) 



Z = Jvifi exp{-5[v5,^*]} 



(2) 



with the action 

1 



2m 



+ 9 E 'P*qi-l3'P*q2+q,^(<l3)'Pqifq2 



Usually, jl is interpreted as the chemical potential /i and 
v{q) as the interaction potential v{q) |3l|9|,[l0|| . 

The intuitive derivation makes use of the coherent 
states introduced by Glauber [ pT|] and can be found in 
Ref. It has been suggested earlier, though, that 

the discretization errors, namely the (AT)^-corrections 
in the exponent, are of importance and cannot be dis- 
carded thoughtlessly |12|. In this work, we present an 



(3) 
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approach which controls these higher order corrections. 

As a consequence, the parameter fl obtains a correc- 
tion from the mathematical manipulations and is to be 
interpreted as iJ + Q, where is a constant depending on 
v{q). The derivation we present here leads to a unique 
relationship between the Hamiltonian and the functional 
integral. 

This issue is related to the work of Christ and Lee 
Q who obtain extra terms in the quantum Yang-Mills 
Hamiltonian from a gauge transformation on the oper- 
ator level. They claim that this is connected to non- 
negligible terms of higher order in the discretization pa- 
rameter e in the functional integral. The existence of 
such possible subtleties in the discretization is our start- 
ing point. We propose a formulation that avoids those 
unpleasant higher order contributions. 

There has been recent significant improvement in the 
understanding of phase space integrals by de Boer, 
Peeters, Skenderis, and van Nieuwenhuizen ^j^, see also 
|l3| . They derive the action and Feynman rules for 
Hamiltonian operators in curved spacetime which are at 
most quadratic in p. This is applied to the non-linear 
sigma model. In this work we are interested in the sta- 
tistical physics applications of functional integrals where 
the interactions typically involve higher powers of p. We 
therefore take a different approach, which allows us to 
treat momentum and location operators on equal foot- 
ing. (There is no restriction to Hamiltonians at most 
quadratic in p). Our simple treatment leads (for flat 
spacetime) to a perhaps unexpected result, namely a cor- 
rection in the action which we can connect to experimen- 
tal findings such as the critical chemical potential at the 
lambda-transition of helium-4. 
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II. FUNCTIONAL INTEGRAL IN QUANTUM 
MECHANICS 



The mathematical issues of path integration are most 
easily studied in one dimensional quantum mechanics, 
which is a finite theory. The generalization of the results 
from this toy model to quantum field theory is straight- 
forward and given below ( IV ) . We consider the canonical 
partition function with the Hamiltonian 



ma) a + —a) a) aa. 



(4) 



We introduce location and momentum operators x — 



+ a) / \/2 and p — 



/\/2, which obey the 



usual commutation relation [x^p] = i. Bringing the in- 
teraction term into a symmetric ordering yields 



H = Ha + Hb — 2" ^ ~8~ 



(5a) 



with 



nA = \ im-X) if+x^) + A {i^ + 2p^'+f) , (5b) 

= i (m-A) if+x') + A (i^ + 2x'p^+p') . (5c) 
4 16 

The term —^X(^x^+p^^ as well as the constant term 
— Y + ^ arise from the use of the commutators in the 
manipulation of (^). The trivial constant term will be 
neglected for notational simplicity. 

This particular ordering is a suitable starting point for 
the formulation of the functional integral. We will show 
that it avoids unpleasant discretization corrections which 
would be present for other formulations like the one based 
on coherent states. 

We denote the location and momentum eigenstates by 
\x) and \p) and follow the standard procedure for writ- 
ing the partition function as a functional integral (with 
e = (3/{2N) and X2n = ^a)'- 



Z = Tr e"'^^ = Z 

N-l 



N 



(6) 

In expanding the exponentials for small e, we will take 
special care of the e^-contributions. The matrix elements 
are evaluated by inserting further location and momen- 
tum eigenstates 

{xk+2\l-2erL + 2e^re\xk) 
dpk+i dxk+i dpk 

X [{xk+2 \pk+i){pk+i\xk+i){xk+i \Pk) {pk\xk) 

-e {{xk+2 \pk+i) {pk+i \'HA\xk+i){xk+i \pk) (pk\xk) 



+ {xk+2 \pk+i){Pk+i \xk+i){xk+i \pk) {Pk ITi-Alxk) 

+ {xk+2 \HB\pk+l) {pk+l\xk+l){xk+l \Pk) {Pk\xk) 
+ {xk+2\Pk+l){Pk+l\xk+l){xk+l\'HB \Pk){Pk\Xk)) 

-f 2e^ {{xk+2 \pk+i){pk+i \HA\xk+i) {xk+i \pk) {pk\T-LA\xk) 

Mxk+2 \'HB\Pk+l){Pk+l \xk+l){Xk+l [HbW) (pk \xk) 
+{xk+2 \pk+l){Pk+l \'HA\xk+l){xk+l iTi-BlPk) {Pk\xk) 
+{xk+2 \HB\pk+l){Pk+l \HA\xk+l){xk+l \Pk) {pk\xk))] ■ 

(7) 

In the above expression, all matrix elements are trivial to 
evaluate, and we see that the e^-term is almost exactly 
the one we need for reconverting the expansion into an 
exponential. All we need to do is adjust some of the in- 
dices (fc -t- 2 ^ A: -|- 1 <-> A:) in this term, but this gives a 
correction of 0{e^) (as will be evident in (^)). Hence we 
can write 



Z 



N 



dxo dx2N-i dpo dp2N-i 



2tt V2tt VStt 

2N~1 



) l^P'' i^k+i-Xk) - I (m-A) {xI+pI) 



k=0 



eA 



{xt+plxl+xl^,pl+pt) +0{e^)\. (8) 



The crucial question is whether Zjv has a well defined 
limit for iV ^ 00 such that the corrections 0{e^) can 
be neglected. Since this is closely related to the issue of 
rapidly varying x and p (as functions of r = /3/c/ (27V)), it 
is useful to perform a Fourier transform. The definition 
of the Matsubara-modes i„, Pn 



N 

Xk — ^ ^ 6 ' ^ ^ Xjij 

n=-N 
N 

n=-N 



^2TTink/{2N)~ 



„27rm(fe+l/2)/(2Af) ~ 



(9a) 
(9b) 



(with i* 



p_„ and xo, Po, xn, Pn real) 



corresponds in the imaginary time language to taking 
the variables p{t) at locations between the x{t). From 
(^), we see that shifting the fc-index introduces a fac- 
tor of 1 -I- 0{1/N). We also turn to the more conve- 
nient language of complex fields by substituting x^ — 
{ip*-„ + Vn) /\/2 and p„ = i (v3l„ - Lp^) / \/2. This yields 



Z 



N 



j Vp exp [-{S + AS) + 0{1/N^)] 



(10) 



with 



2 



5 = ^ [2TTin(plLpn + /? (m - A) (p* ( 



f3X 



\ " * * 



(11) 



A5 = ^[4.A.sin(^)-2. 
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E 



1 — cos 



7r(ni + 712) 
2iV 



x('^-^2^ - '/'n2+n)('^-^2 ~ '^«2)('y'-^u - "Pn^)- (12) 

We notice that the couphng constants in the A5- 
correction are cx ©(l/TV^). (The n-dependence is not 
important in this respect. Up to a constant prefactor of 
Z, the contributions of the high n Matsubara-modes are 
effectively cut off by suppression factors oc from the 
squared propagator, which arise from the first term in 
S.) Eq. (|l^) is therefore particularly suited for an un- 
equivocal definition of the functional integral as the limit 

— !■ 00 is taken. We extend the n-summation from —oo 
to 00 in the action S and drop A5 as well as the other 
O ( 1 /iV^ )-corrections . 

In the "imaginary time language" (using yj(T) = 
g2irmT//3j^^-j^ the action corresponds to 



S- 



li/2 



dr 



-13/2 



(13) 



In this simple case, the only difference to the naive use of 
coherent states is the shift of A in the mass term (besides 
constant terms). 

Only the symmetric ordering of (||) avoids unpleasant 
1 /A^-corrections. This criterion leads to a unique trans- 
lation into a continuous functional integral. We can also 
take the perspective from the problem of operator order- 
ing [p^|~^. To think of the easiest example, the different 
operators xp and px give the same zeroth order contribu- 
tion in the action. This leads to an ambiguous functional 
integral if higher order corrections are not taken into ac- 
count. The difference between these operators must oc- 
cur elsewhere in the functional integral, and we see that 
it precisely appears in the l/A^-terms. These can give 
finite contributions if the evaluation of the functional in- 
tegral preceeds the limit iV — > cx). In the transition to 
the continuous functional integral, they would thus be 
erroneously neglected. 



III. A SIMPLE TEST 

In order to check our assertion that the quadratic term 
should be TO — A instead of to, we calculate the ther- 
mal equilibrium value of the momentum squared from 
the functional integral in first order perturbation theory. 
This is compared to a direct evaluation from ordinary 



quantum mechanics. The two conclusions we draw from 
this will be that the functional integral gives finite results 
and therefore needs no regularization and that only the 
quadratic term to — A gives correct quantitative results. 

We start with the simple quantum mechanical calcu- 
lation. We denote the eigenstates of the unperturbed 
Hamiltonian operator (A = 0) by \l). With {l\p'^\l) = 
I + ^, we obtain the thermal expectation value 



s2\ 



-fJml-^f3\l{l-l) 



-I3ml-^I3\1{1-1) 



(14) 



Due to the fast convergence of the two sums, this can 
easily be evaluated numerically. 

In the functional integral, the expectation values of 
operators are easily computed by adding to TL source 
terms like Kp^, e.g. (p^) = d\nZ{K)/dK\j^_^. As long 
as the source terms do not involve products of p and x, 
they simply add to the action the corresponding terms 
with the replacements f{p) J dr f (i {(p* — ip) /V2) 
and f{x) J dr / (((^* + 1^9) /\/2) . One obtains 

(p^) — ~5((v* ~ f)'^) with the usual definitions of ex- 
pectation values in the functional integral {0{(f,p*)) ~ 
Z^^ J Dp 0{(p, (p*) exp{—S). A standard calculation in 
first order perturbation theory yields 



(p^) = - coth 

' 2 2 



/3Acoth^ 
4sinh2^ ' 



(15) 



where M — m corresponds to the naive coherent state 
approach and M = to — A is our suggestion. 

In table ||, we show how this compares for both choices 
of M and different values of to and P with the exact re- 
sult. We take A to be small against T = 1//3 and to 
in order to justify the perturbative approach. The re- 
sults clearly demonstrate that the shift in the mass term 
is necessary for quantitatively correct results, and this 
alone should be convincing that the correction we intro- 
duced to the standard functional integral is necessary for 
the correct interpretation of the parameters in the action. 

TABLE L Momentum squared (p^) for two different sets 
of f3 and m by exact quantum mechanics, eq. (|l^), and by 
first order perturbation theory from the functional integral, 
eq. (p^. We compare the suggested shift in the mass term 
(M = m— A) with the naive coherent state approach (M = m). 
All entries must be divided by 10^, and the zero-temperature 
value of 0.5 is subtracted from (p^). 



A 


for /3 = l,m 
exact Al — m — X 


= 5 
M^m 


(p^) for /3 = 3, m 
exact M — rn — \ 


= 1 
M = m 





6.7837 


6.7837 


6.7837 


52.396 


52.396 


52.396 


2 


6.7835 


6.7835 


6.7698 


52.361 


52.359 


52.030 


4 


6.7833 


6.7832 


6.7560 


52.327 


52.320 


51.665 


6 


6.7831 


6.7830 


6.7421 


52.293 


52.278 


51.299 


8 


6.7829 


6.7827 


6.7283 


52.259 


52.232 


50.934 
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IV. FUNCTIONAL INTEGRAL FOR 
INTERACTING BOSONS 

The preceeding derivation for quantum mechanics can 
directly be generahzed to an interacting ensemble of 
bosons and therefore becomes relevant to the treatment 
of the superfluid transition of hehum-4 or Bose con- 
densation. We describe the system through the Hamil- 
tonian operator (|l|), assuming the existence of a UV- 
momentum-cutoff A, < A^. The generahzation of our 
procedure gives the action (|3|) with (apart from constant 
terms) 'O(g) = v{q), p, = + fl, and 



n 



2 ^ 



w(0) + v{q) 



(16) 



We obtain a cutoff-dependent "counterterm" to the 
chemical potential, which will be of importance for com- 
parison with experiment. As it turns out, this cancels 
a similar term generated by the one loop approximation 
for the "full" inverse propagator. This cancellation leads 
to the correct low-temperature phononic dispersion rela- 
tion. 

We have calculated the chemical potential at the A- 
transition of helium-4 under vapor pressure conditions 
by means of a truncation of an exact renormalization 
group equation iQ. The complete results are presented 
elsewhere jl9|. Here we just refer to the value of the criti- 
cal chemical potential, which comes out to be —6.75 K in 
good agreement with the experimental value of around 
—7.4 K |20|. As il in our case has a value of approxi- 
mately 12 K (with a cutoff given by the atomic length 
scale), neglecting this correction would have resulted in 
a strong deviation from experimental values. 



V. CONCLUSION 

In this letter, we advocate a specific approach to 
the derivation of the functional integral for an interact- 
ing quantum model and its generalization to interacting 
bosonic systems. We are led to a shift in the mass term of 
the action as compared to conventional approaches based 
on the use of coherent states. A perturbative calculation 
of the expectation value of the squared momentum in the 
interacting quantum system and a renormalization group 



calculation of the critical chemical potential for the A- 
transition of helium-4 confirm the necessity of this shift 
and strengthen our faith in the suggested prescription. 

These findings should be important to renormalization 
group treatments of statistical systems such as Bose- 
Einstein-condensation for interacting systems as well as 
to the understanding of functional integrals in general. 

We are grateful to Andrew Waldron for drawing our 
attention to the works of Christ and Lee and on the treat- 
ment of the non-linear sigma model. 
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